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ABSTRAC  r  (ConliiHM  an  rAYAfAA  AIRA  II  AAAtAAAAr  anR  IRaiKIIA  »a  *IacA  imaaSaa) 

The  approximate  expression  for  the  growth  rate  of  the  transverse  electric  tTEI  minle  with 
the  azimuthal  mode  number,/,  the  radial'.node  number  n,  and  the  magnetic  harmonic  niimbei  t 
of  the  electron  cyclotron  mater  instability  with  zero-tcmpertured  beam  electrons  is  dciived  and 
analytically  examined.  It  has  been  found  that  there  exist  two  branches  of  instabilities  drstinguisbed 
by  their  axial  wavenumbers;  the  short  and  the  long  wavelength  modes,  and  the  estimates  of  their 
growth  rates  have  been  obtained.  With  maximization  processes  of  the  growth  rates,  the  optimum 
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20.  Abstract  (continued) 
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'conditions  on  various  parameters  have  been  found, 
with  the  beam  center  coinciding  with  the  axis  yields 
behavior  of  the  growth  with  respect  to  the  magnetic 


In  particular,  it  has  been  found  that  s  =  I  mode 
the  highest  growth  rate.  The  exponential 
harmonic  number  s  is  also  obtained. , 
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ANALYTIC  ANALYSIS  OF  UNEAR  GROWTH  RATES  OF  GYROTRON  FOR 
TE|„s  MODE  WITH  ZERO-TEMPERATURE  BEAM  ELECTRONS 

I.  INTRODUCTION 

The  general  dispersion  equation  for  transverse  electric  (TE)  modes  of  the  electron  cyclo¬ 
tron  maser  instability  was  previously  obtained  in  Ref.  I  (hereinafter  referred  to  paper  A).  It  is 
the  objective  of  this  paper  to  analytically  examine  the  dispersion  relation. 

In  the  spirit  of  application  of  the  cyclotron  maser  instability  to  gyrotron  amplifier  devices, 
we  will  identify  different  instability  types  and  will  find  optimum  conditions  in  experimental 
parameters  that  maximize  the  growth  rates.  The  analytical  expression  for  the  growth  rales  is 
obtained  with  the  tenuous  beam  approximation  that  is  consistent  with  the  assumption  in  deriv¬ 
ing  the  general  dispersion  equation  (see  Chapter  II  in  paper  A). 

Although  the  dispersion  Eq.  ((70)  in  paper  A)  is  valid  for  arbitrary  distribution  function 
of  beam  electrons,  we  will  restrict  ourselves  to  zero-temperature  beam,  whose  dispersion  rela¬ 
tion  is  given  by  Eq.  76  in  paper  A.  Also  we  will  consider  the  gyrotron  with  single  conductor, 
for  simplicity,  even  though  some  of  the  results  are  still  valid  for  coaxial  gyrotron  devices. 

The  geometry  of  the  gyrotron  we  consider  is  that  of  cylindrical  wall  within  which  the  elec¬ 
tron  beam  is  passed.  The  cross  section  of  the  device  is  shown  in  Eig.  1.  The  conducting  wall  is 
located  at  z  —  Af,  and  the  beam  electrons  of  constant  Larmor  radii,  r,  are  gyrating  with  constant 
guiding  center  radii,  z„  around  constant  magnetic  Reid  The  axial  motion  of  electrons  are 
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characterized  by  constant  velocity  v_..  The  distribution  function  describing  such  electrons  is  a 
6-function  given  by  Eq.  (73)  in  paper  A,  and  the  resulting  dispersion  relation  for  TE,„,  mode 
becomes  Eq.  (76)  in  paper  A.  Here  mode  identification  integers  /.  n.  and  s  represent  the  mode 
numbers  for  the  azimuthal  (W),  the  radial  (r)  directions,  and  the  magnetic  harmonic  number  of 
Doppler-shifted  cyclotron  motion,  respectively.  In  addition,  the  axial  mode  is  identified  by 
wavenumber  k. 

The  dispersion  relation  eq.  (76)  in  paper  A  is  obtained  under  the  tenuous  beam  approxi¬ 
mation  within  the  framework  of  linearized  Vlasov-Maxwell  equation  with  form  of 
explikz  -(■  i/W  -  luii].  The  tenuous  beam  assumptions  can  be  written  as 

I'/y,)  <  <  <  1 .  ( 1 ) 

where  the  Budker  parameter  r  is  defined  by 

■Vc^ 

~  .  (2) 

and  y,,  is  the  total  relativistic  factor  of  the  electron.  Here  the  total  number  of  electrons  per  unit 
axial  length  /V  is  given  by 

;V  =  27r  J"  nir  /„(y,  p..P„).  (3) 

The  procedures  of  deriving  the  dispersion  relation  is  very  complicated,  and  we  refer 
details  to  paper  A.  It  is  sufficient  to  say  that  the  tenuous  beam  assumption  (I)  allows  us  to 

neglect  self-field  effects,  and  that  usage  of  constants  of  electron  motion  (y.  />.,  /*«)  as  the  coor¬ 

dinates  for  beam  distribution  function  /,i  eliminates  much  of  the  mathematical  complexity. 
Here  the  constants  of  motion  y,  />.  and  refer  to  the  total  relativistic  factor  of  the  electron, 
the  linear  momentum  in  axial  direction,  and  the  azimuthal  canonical  angular  momentum. 

After  brief  discussions  of  the  dispersion  equations  in  Chapter  If,  we  identify  two  different 
types  of  instability  (short-and  long-wavelength  modes),  and  derive  approximate  expressions  for 
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their  growth  rates  (Chapter  II).  From  Chapter  IV  and  on,  we  limit  our  attention  to  more 
interesting  long- wavelength  modes.  The  stability  boundary  is  found  in  Section  A  of  Chapter 
IV,  and  the  estimates  of  growth  rates  and  the  matching  conditions  on  wall  location  are  obtained 
in  Section  B.  Further  optimization  of  the  growth  rates  with  respect  to  various  parameters  are 
given  in  Chapt.  V. 


It  should  be  emphasized  here  that  our  intention  in  this  paper  is  to  obtain  general  tendency 
of  the  dispersion  equation  without  help  of  numerical  calculations.  For  this  purpose  we  fre¬ 
quently  resolve  the  graphic  solutions. 


11.  DISPERSION  EQUATION 


The  dispersion  relation  for  TE,„  mode  for  the  device  described  in  Chapter  1  is  given  by 


Eq.  (76)  of  paper  A.  Namely 


where 


The  various  terms  in  Eq.  (S)  are  dellned  by 
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and  (he  functions  H  =  H,i{Xq.Xi),  Q  =  Q^iiXo,  Xl),  and  U  =  Xi,)  are  given  by 

H^(Xo.Xl)  -  lJ,.,iXo)JAXL)]\ 


QAXo.Xl)  -  J,.,(X^  JAXi,) 


s{s^-X^) 


J,.,{Xc)J,(Xl) 


XHs  -  /)  .  .  .  , 

+  - T.  J>(Xi) 


Xo 

LI„(Xo.Xl)  -  J,(Xl)  J[(X,,) 


(2s^  -  Xl)  XiM-lV  ,  .j 


Xl 


Xl 


JUXo) 


-  JTi  U-,(Afo)  . 


The  arguments  of  H.  Q,  U  are  defined  by 


where 


„  _  x,„ 

^0  =  '’o- 

X,  =  rt. 


Ji(Xi„)-0. 


(7) 


(8) 


(9) 


The  solution  of  Eq.  (4).  that  is,  if  D(w  -  <i>,  +  /«,)  -  0,  determines  the  eigen  frequency 
w,  and  the  growth  rate  |(i>,|  if  oi,  ^  0.  The  first  term  in  Eq.  (S)  describes  free  wave  guide 
motion  with  eigen  frequency  oti„,  while  the  second  terms  represents  the  correction  to  free  wave 
guide  equation  due  to  the  cyclotron  motion  (with  frequency  O,)  of  beam  electrons.  It  is  thus 
the  second  term  that  gives  rise  to  the  instability  with  m,  ^  0. 


For  systematic  analysis,  we  normalize  quantities  as  follows 


Ui 

n. 

-  s«,./yo 

W/n 

A, 

R  - 


cX,„ 


y,su,/y„ 
'•o  -  ^  '0- 


(10) 
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where  the  axial  relativistic  factor  y,  is  defined  by 


y,  5  (1  ^ 


With  thus  normalized  guantities,  the  dispersion  equation  D  —  0  can  be  written  as 


DiH:  k) 


'fix' HO, 


(w  -  0,) 


where 


n.  =  kfi,  +  1. 

<U) 

Xl  =  Pi,  A’o  sXi„ff) 

Here  Ai„  is  small  quantity,  typically  in  the  order  of  10"’. 

The  normalized  dispersion  function  D(i»,k)  is  convenient  in  that  the  magnetic  harmonic 
number  s  is  factored  out  except  for  the  complicated  functions,  H,  Q,  and  U.  Hereinafter,  we 
will  use  the  normalized  quantities  (Eqs.  (12)  and  (13))  with  tilde  dropped.  When  we  consider 
the  real  parameters,  the  normalization  Eq.  (10)  will  be  referred  to. 


The  waveguide  eigenfrequency  «/,,  and  the  Doppler-shifted  cyclotron  frequency  ft,  as 
functions  of  axial  wave  number  are  illustrated  in  Pig.  2.  It  is  useful  to  define  the  following 


wavenumbers; 


n,  -  -  0  at  k  -  k*4. 
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n*  -uii-O  at  k  ~  k%. 
and 

n.  -  0  at  A  -  A,,. 

It  is  easy  to  show  that,  from  Eq.  (13), 

A"*  -  y/{^,  ±  11 

A-,  -  y;(/8,  ±  (1  -!//?•) 
and 


(15) 

(16) 

(17) 

(18) 


(19) 


The  waveguide  fretiucncy  a>/,  is  a  hyperbolic  function  of  A  with  the  cutoff  frequency 


R 


(i.e. 


intersection  of  cu,,  with  lu-axis),  and  the  cyclotron  frequency  fl,  is  a  straight  line  of  slope  /3. 


intersecting  <o-axis  at  1.  U  is  evident  from  Eq.  (18)  that  for  /(  <  1,  wi„  is  always  greater  than 


n,,  and  particularly  for  /?  -  1,  fl,  is  tangential  to  <u/„.  Figure  2  shows  an  example  for  R  >  1. 


Before  we  obtain  the  expression  for  the  growth  rate  w,  from  Eq.  (12),  we  will  first 
demonstrate  the  presence  of  the  complex  roots  (i.e.  to,  ^  0)  in  Eq.  (12).  We  rewrite  the 
dispersion  equation  as 


where 


D  “  D(i  +  f)|  +  Dy  +  27 j. 


(20) 


3 

1 

3 

III 

Q 

(21) 

(22) 

^  R^  («-n,)  • 

(23) 

=  -  -  film. 

(24) 
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Here  small  guantity  i)  is  defined  by 


n  =  — ^  HAXM. 

and  the  typical  magnitude  of  17  is  in  the  order  of  10  *. 


(25) 


For  given  parameters  and  k  the  dispersion  function  £)  is  a  function  of  <»,  and  the  roots  of 
D(o»)  -  0  will  determine  the  eigenfreguency.  In  order  to  have  growing  mode,  the  eguation 
D(ot)  -  0  must  have  complex  roots.  We  note  that  in  the  absence  of  beam  electron  (i.e. 

-  0  -  »j),  the  dispersion  eguation  0-0  becomes  Do  -  0.  The  free  waveguide  function 
Do(w)  is  a  quadratic  function  of  w,  and  the  free  waveguide  eguation  Do  —  0  has  roots  at 
Cl*  -  ±  o*(,,  representing  freely  (without  growth  or  attentuation,  i.e.  <0,  -  0)  travelling 
waveguides  in  opposite  directions.  With  electrons  (v  ^  0).  we  can  easily  show  that,  due  to  the 
resonance  term  (&*  -  fl,)  in  denominator  of  D|  and  Oj,  and  the  smallness  of  the  coefficients 
of  0|,  Z>2,  and  Oj,  the  dispersion  function  D(<i*)  can  be  approximated  as 


at  jci*  —  ft,  I  >  >  0,  and 


D(u>)  ==  Dpicu) 


(26) 


D(w)  —  Do(ci*)  +  /)|(<u)  (27) 

near  &»  -  ft,,  where  the  magnitudes  of  Dj  and  O3  are  much  less  than  that  of  £>1  (see  Egs.  (22) 

-  (24)).  Since  Do(ut)  -0  (Eg.  26)  have  real  roots,  and  complex  roots  are  possible  only  for 

Do(iv)  +  Di(io)  —  0  (Eg.  27),  we  conclude  that  instability  may  be  present  only  near  w  —  ft,. 

We  therefore  examine  the  behavior  of  Doito)  +  D,(a»)  with  respect  to  <1*  near  ft,,  which  are 

illustrated  in  Figs.  (3a)  and  (3b).  At  two  representative  values  of  k  (noted  as  I  and  II  in  Fig. 

2),  their  approximate  dispersion  functions  Dofeu)  +  D|(fi>)  with  Do(u>)  for  comparison  purpose 

are  drawn  near  a*  -  ft,  in  Fig.  3(a)  and  (b).  Note  that  k'  <  k1  <  k",  and  that  at  k  -  k' 

(Fig.  3a)  ft/  >  Cl*,',  while  at  k  -  k"  (Fig.  3b)  ft,"  <  a*,"  It  is  apparent  from  Fig.  3  (a)  and 

(b)  that  the  difference  of  dispersion  function  with  electron  beams  (Do  +  0|)  from  free  wave 


CHOE  AND  AHN 


guide  function  (D,,)  is  localized  near  oi  -  n,.  The  combined  dispersion  function  hardly 
changes  the  free  waveguide  eigenfreguency  cu,,,  while  at  the  opposite  side  of  Jl,  from  it 

IBD  I 

wj  -r—  -  0  and  D^—r  >  0  .  As  can  be  seen  in  Appendix  A,  the 

oot  ) 

fact  that  the  function  D(w)  has  a  local  extremum  near  <i>  -  H,  means  that  the  dispersion  equa¬ 
tion  Z)(cij)  -  0  yields  complex  solutions  (a>,  ^  0)  near  that  point.  Since  imaginary  part  of  the 
solution  of  Diw)  -  0  represents  the  growth  rate,  the  presence  of  the  local  extremum 
corresponds  to  instability  near  oj  -  11,.  The  presence  of  the  local  extremum  near  u  -  11,  is 
due  to  the  fact  that  D|  is  proportional  to  (<o  -  11,)'%  which  yields  steep  slope  near  oj  -  11, 
and  same  signs  as  it  crosses  lu  -  11,  line.  We  have  graphically  shown  that  the  dispersion  equa¬ 
tion  D(<o)  -  0  can  have  complex  solutions  and  thereby  growing  modes  near  oj  -  11,.  In  the 
next  chapter  we  will  hnd  the  approximate  expression  for  the  actual  growth  rate. 

III.  GROWTH  RATE  EQUATION 

The  dispersion  Eq.  (12)  can  be  written  as 
A-  1(11,’  -  OJ,;)  -t-  211,  A  +  A^i 

+  /8i^//(n,2-A2)  -2  a  (29) 

A  5  1 

-  AMU  -  /3/  H)  -  0. 
where  the  small  frequency  shift  A  is  defined  by 

A  =  o)  -  n,  =  A,  -t-  /  A,.  (30) 

When  the  tenuous  beam  is  assumed  (Eq.  (D),  it  can  be  shown  later  that 

|A|  «  1.  (31) 

We  can  find  the  analytical  solutions  of  Eq.  (29)  for  different  regions  of  k  where  some  approxi¬ 
mations  can  be  made  on  the  terms  in  Eq.  (29) 
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A.  Short  Wavelength  Mode  (k  <  kt.  k  >  *+) 

For  k  <  kt  OT  k  >  ki  (see  Fig.  2  and  Eqs.  (14)  and  (17)),  we  have 

»  2n,  iA|. 

»  lA'l. 

H \{l^  -  k^  »  2  ^  HQ,,  (32) 

\iU-k^  »  lAlMt/ 

SO  that  Eq.  (29)  can  be  approximated  or 

AMn;  -  C/5,)  + (n.^  -  A-^) -0.  (33) 

H 

Here  small  quantity  17  is  defined  in  Eq.  (2S).  The  solution  of  the  quadratic  Eq.  (33)  is  given  by 


A,  -  0, 


where  A,  and  A,  are  the  real  and  imaginary  parts  of  A,  and  c/„  and  fl,  are  given  by  Eq.  (13). 

I 

We  note  that  the  growth  rate  A,  is  proportional  to  (17)  that  it  is  a  monotonic  function  of  axial 

I 

wavenumber  k,  approaching  ri^/R  as  |A'|  increases  to  infinity,  and  that  the  growth  rate  is  a 
monotonically  decreasing  function  of  the  wall  location  R.  These  are  illustrated  in  Fig.  4.  it  can 
be  easily  proved  that  al  k  —  k±  (see  Eqs.  (14)  and  (17)),  there  exists  no  instability  since  the 
dispersion  equation  becomes  linear. 

B.  Long  Wavelength  Mode  {kt  >  k  >  kt)- 
For  kt  >  k  >  kt,  (see  Fig.  2,  and  Eqs.  (14)  and  (17)),  we  have 

n,  »  |A|. 

P^^H  (n?  -  k^)  »  2|Al  -2  -  13/  wn, ,  (36) 

5 

P^H  (n/  -  k^)  >S  |A^|  \  U  -  p^^H\. 
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SO  that  tig.  (29)  can  be  approximated  as 


A'  +  (i  -  0. 

Here  guaniilies  a  and  h  are  defined  by 


a  = - 


(37) 


t>  = 


2H. 

=  JL 


(38) 


(39) 


R-  211. 

Note  that  we  deliberately  retain  A'  term,  since  the  long  wavelength  mode  (A'’  >  k  >  k^) 
included  the  A-region  where  the  guantity  a  vanishes  (Eg.  (IS)  and  Eig.  2).  The  cubic  Eg.  (37) 
can  be  formally  solved  as  shown  in  Appendix  B.  Namely,  the  complex  solution  of  Eg.  (37)  is 
given  by 


A,  -  y  +  /fi  I. 
fx  I  I 

A,  -  ^  I. 


where 


~  2  27  * 


(l‘  =  h 


4  27 


The  necessary  and  sufficient  condition  for  instability  (i.e..  A,  >  0)  is  readily  obtained. 

(/'  >  0. 

We  note  that  for  the  long  wavelength  mode  (A!1  >  A  >  A^), 

^  >  0. 


(40) 

(41) 

(42) 

(43) 

(44) 


H,  >  B.>  0. 
B,B  - 


(45) 


The  growth  rate  of  the  long  wavelength  mode.  Eg.  (41)  will  be  extensively  analysed  in  the  next 
two  chapters. 
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We  have  seen  that  there  exist  two  branches  of  the  electron  cyclotron  instabilities.  These 
are  identitieU  as  the  long  wavelength  and  the  short  wavelength  modes.  These  modes  are 
separated  by  band  of  stable  region  in  the  axial  wavenumber  (A)  near  A  —  A'J  (see  Fig.  4>.  and 

I 

the  growth  rate  of  the  long  wavelength  mode,  as  it  turned  out  later,  is  in  the  order  of  rj '  (see 

En.  (25)),  while  that  of  the  short  wavelength  mode  is  in  the  order  of  .  At  usual  parameters 
(see  Eq.  (50)),  the  growth  rate  of  the  long  wavelength  mode  is  four  times  higher  than  that  of 
the  short  wavelength  mode.  We  also  note  that  the  axial  group  velocity  of  the  waveguide  for  the 
long  wavelength  mode  is  in  the  order  of  the  beam  axial  velocity,  while  that  for  the  short 
wavelength  is  much  higher  and  is  in  the  order  of  speed  of  light  (<•).  The  combined  elVects  arc 
such  that  the  ampliftcation  of  the  long  wavelength  mode  is  much  more  dominant  than  that  of 
the  short  wavelength  mode  when  both  are  present.  It  should  also  be  emphasized  that,  as  dis¬ 
cussed  in  the  next  chapter,  there  exist  parametric  optimum  conditions  that  maximize  the 
growth  rate  of  the  long  wavelength  mode,  while  such  conditions  does  not  exist  for  the  short 
wavelength  mode,  which  makes  it  difticult  to  design  an  elTtcient  experiment  utilizing  the  short 
wavelength  mode  It  is,  however,  possible  that  for  some  particular  experimental  situations  the 
amplification  due  to  the  short  wavelength  mode  may  be  comparable  to  or  even  higher  than  that 
due  to  the  long  wavelength  mode.  For  this  reason,  we  can  not  completely  forsake  the  study  on 
the  short  wavelength  mode.  It  is  pointed  out  that  there  is  neither  the  experimental  proof  of  the 
existence  nor  the  plausible  explanation  of  the  physical  mechanism  responsible  for  the  short 
wavelength  mode  yet. 

IV.  LONG  WAVELENGTH  MODES 

The  long  wavelength  mode  instability  is  detined  as  a  growing  mode  whose  axi-j 
wavenumber  is  within  the  range  of  A''  <  A  <  A^,  where  arc  detined  by  Eq.  (17).  The 
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growth  rate  for  the  long  wavelength  mode  is  given  by  approximate  solution  (41).  Generally  the 
growth  rate  A,  is  a  function  of  following  normalized  parameters. 

A,  -  A,  (A,  r.,;  n.  I.  s;  ft,  fi,.  v)  (46) 

The  various  parameters  are;  the  axial  wavenumber,  the  conducting  wall  radius,  the  beam  center 

radius,  the  radial  node  number,  the  azimuthal  mode  number,  the  magnetic  harmonic  number, 
the  transverse  speed,  the  axial  speed,  and  the  density  (Budker  parameter).  In  this  chapter  we 
will  discuss  the  behavior  of  the  growth  rates  with  respect  to  the  axial  wavenumber  (A)  and  the 
wall  radius  (R),  while  all  the  rest  parameters  are  assumed  to  be  given  The  examinations  con¬ 
sist  of  analysis  on  properties  of  quantities  appearing  in  growth  rate  Eq.  (41),  the  determination 
of  the  instability  boundary,  and  maximization  of  the  growth  rate  in  k  -  R  space.  We  adopt  the 
notation  "  as  optimum  value  that  maximize  the  growth  rate. 


It  is  our  intention  to  obtain  as  much  information  on  the  general  properties  and  tendencies 
of  the  growth  rate  as  possible  without  actually  solving  the  growth  rate  Eq.  (41).  This  analytical 
approach  is  possible  primarily  due  to  the  tenuous  assumption  (see  Hq.  (I)).  This  small  density 
assumption  can  be  expressed  as 

T)  «  1,  (47) 

where  17  is  defined  by  Eq.  (25),  namely 

2v  V?  lii^ 

T,=  - IJ,.,{X„)  J:{X,)]\  (48) 

yoNi„ 

Here  (Af„,  X/  )  of  Eq.  (25)  is  explicitely  shown,  and  quantities  A',„  and  X,  are  given  by 


■^0  —  '’()• 


Using  typical  experimental  parameters  given  by 

p  -  3  X  10-’.  ft  -  0.3.  li.  -  0.2. 
.V  —  /  —  «  —  1 ,  r„  -  0,  A?  —  1 . 


(8) 

(13) 


(49) 
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the  values  of  small  quantity  r)  and  its  fractional  powers  are 

I 

t)  -  2.7  X  10  ^  7,^  -  •  <*  X  (50) 

I  I 

T,  '  -  6.4  X  10'^  rj"  -  2.5  x  10  '. 

It  should  be  mentioned  here  that  the  mode  s  -  /  -  n  -  1  with  r„  -  0,  R  =  I  yields  the 
highest  growth  rate  in  entire  (/,  n.  s)  modes.  That  is  to  say  that,  for  given  »>,  fii  and  /8,.  the 
values  given  in  Eq.  (50)  are  their  upper  limits.  The  smallness  of  tf  enables  us  to  approximate 
equations  so  that  their  graphic  solutions  and  estimations  of  the  numerical  values  of  solutions 
are  frequently  possible  to  obtain. 

It  is  convenient  to  consider  the  growth  rate  A,  as  a  function  of  the  quantities  a  (A.  ft)  and 
MA.  ft)  dctined  by  Eq.  (38)  and  (39)  respectively.  The  quantity  a  determines  the  difference 
between  two  eigenfruquencies  of  the  system,  namely,  the  Doppler-shifted  cyclotron  frequency, 
1),  (Eq.  (13))  and  the  free  wave  guide  frequency,  iu,„  (Eq.  (13)1;  while  the  quantity  ^  meas¬ 
ures  the  difference  in  phase  velocities  between  the  beam -waveguide  eigenmodes  and  the  free 
electromagnetic  wave  in  vacuum.  We  will  examine  the  functional  dependencies  of  a  and  ^on  A 
and  ft  for  future  references.  With  some  algebra,  the  quantities  a  and  hcan  be  written  as 

a - L.  _  A)  (*  _  A«).  (51) 

2yl  llj 

6  -  -J-  -1-  —  ikl  -  A)  (A  -  kt),  (52) 

2>/  n, 

where  A'i  and  A^  are  given  by  Eqs.  (18)  and  (17).  If  we  differentiate  a  and  6  with  respect  to 
A,  it  is  straightforward  to  show 

(A;,'  -  A)  (A  -I-  A',').  (53) 

dA  2y?  O’ 

-  J-i-  -3-  -L.  (A^-  A)  (A  +  AC).  (54) 

dA  2y?  ” 
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where 


jjjl  =  i  Ir,  (I  +  r-  «?/«’)"'  Tl).  (55) 

hi  ^  i  ^  "•  "6) 

The  plots  of  A  4,  A",  A^,  and  A,J  in  k  ~  R  space  are  given  in  Fig.  5.  We  note  that  A®  >  A*  and 

that  ki  <  ki  K  k%  ^  k1  when  >  1.  For  /?  <  1,  A®±  are  not  defined.  While  A®^  and  A® 

are  functions  of  R,  the  quantities  A^  and  A, fare  independent  of  R.  The  Eqs.  (51)  to  (54)  indi¬ 
cate  that 


and 


a>  0  when 

< 


6  >  0  for  ki 


k%  >  k  >  kt  1 
A  >  Ai  or  A  <  A® I 
>  k  >  kt. 


da 

dk 


>  0  when  A 

< 


(57) 

(58) 

(59) 


>  0  when  A  <  A,f.  (60) 

dk  <  > 

Note  that,  from  Eq.  (58),  6  >  0  for  long  wavelength  modes.  For  future  reference  we  define 
the  wall  radius  of  the  intersection  of  Af  line  and  A"  curve  as  R„  as  shown  in  Fig.  5,  and  its 
magnitude  is  given  by 


With  /a,  -  0.2,  R„ 
(Appendix  C) 


Rn  - 


r.-  +  l 


1/2 


(61) 


1.005.  The  wall  radius  (R)  dependences  of  a  and  b  can  be  shown  as 


da  _  2  1 

dR  ~  n,R^’ 


db 

dR 


-  V 


(n/-A^)  J,(Xl) 


(62) 


(63) 


n,/?^  X,,  j,'(x,,) 

The  determination  of  the  sign  of  db/dR  is  complicated  since  the  sign  depends  on  the  signs  of 
the  terms,  (s^  -  2(/),  J^iX^)  and  J,'  (X^).  It  can  be  shown,  however,  that 
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—  <  0  When  R  >  y.  ((,4) 

as  shown  in  Appendix  C.  H  should  be  noted  that  the  condition  R  >  y.  is  easily  satisfied  for 
weakly  relativistic  (y,,  —  1)  beam  electrons  as  in  usual  gyrotron  devices.  The  derivatives  of  the 
quantities  a  and  h  with  respect  to  the  wavenumber  k  and  the  wall  radius  R  (i.e.,  Egs.  (53), 
(54).  (62),  and  (63))  are  computed  for  the  purpose  of  finding  the  location  of  the  maximum 
growth  rate  in  k  -  R  space  later. 

A.  Instability  Boundary 

Before  discussing  the  magnitudes  of  the  growth  rates,  it  is  instructive  to  examine  the  ins¬ 
tability  criteria.  This  will  determine  the  instability  boundary  in  k  -  R  space,  and  thereby  the 
instability  area,  which  is  closely  related  to  the  bandwidth  and  the  cut-off  value  of  the  wall 
radius,  which  is  important  in  designing  gyrotron  devices. 

The  necessary  and  sufficient  condition  for  the  long  wavelength  instability  is  given  by  Eg. 
(44)  (see  Appendix  B).  which  can  be  written  as 

</  '  ^  b 

27  >  4  •  (65) 

The  quantities  a  and  A  are  given  by  Egs.  (51)  and  (52),  We  again  emphasize  that  6  is  positive 

but  small  (proportional  to  fof  the  long  wavelength  modes.  We  realize  that  solving  Eg.  (65) 
is  equivalent  to  solving  the  growth  rate  Eg.  (41).  and  resolve  to  graphic  solution  of  Eg.  (65) 
avoiding  much  of  algebra.  The  solution  of  the  stability  boundary  Eg.  (65)  is  graphically  illus¬ 
trated  in  Fig.  6.  The  plots  of  (<iV27)  and  (-  A/4)  versus  the  wavenumber  k  arc  given  for 
several  values  of  the  wall  radii  R  designated  as  I  to  V.  The  numerical  values  of  R  decreases  as 
the  notation  changes  from  1  to  V.  such  that  R,  >  R„  >  R,„  -  1  >  /?„  >  /?,.  The  instabil¬ 
ity  region  in  k  is  bordered  by  two  intersecting  points  of  (-  A/4)  and  (<i  V27)  curves  for  each  R. 

Since  the  curve  of  (fl  V27)  at  given  R  nests  all  (<j  V27)  curves  at  smaller  wall  radii  (see  Eg. 


15 


CHOE  AND  AHN 


(51) )  inside,  and  (—  b/A)  is  a  weak  function  of  R  (the  coefficient  is  small  number  t),  see  Eg. 

(52) ),  we  conclude  that  the  instability  range  AKC/t)  is  a  montonically  increasing  function  of  R, 
namely 

Ik  («,)  >  IkiR})  if  «,  >  «2,  (66) 

and  that  there  exists  a  cutoff  wall  radius  Rg  such  that 

XiR  <  Rg)  -  0.  (67) 

These  features  are  apparent  from  Fig.  6.  The  distance  between  two  intersecting  points  (i.e., 

AA)  becomes  smaller  as  I  goes  to  V  (i.e.,  as  R  decreases),  and  for  small  R  (e.g.,  R\/)  the  two 

curves  do  not  intersect,  corresponding  no  instability.  The  schematic  graph  of  the  stability 

boundary,  denoted  as  kg,  \n  k  —  R  space  is  drawn  in  Fig.  7.  The  instability  occurs  only  to  the 

right  of  the  boundary  curve  kg,  and  we  note  that  kg  intersects  with  k“  at  Rg,  the  cutoff  wall 

radius,  since  ia^/27)  is  maximum  at  k,‘l  (see  Eg.  (S3)). 
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and 

1  -  y  (70) 

With  parameters  given  by  Eqs.  (49)  and  (50),  the  cutoff  wall  radius  Rb  is  found  to  be  ~0.904. 

The  bandwidth  of  the  gyrotron,  if  we  define  it  by  the  range  of  the  axial  wavenumber  k 
where  the  instability  exists,  is  given  by  Eq.  (69).  With  this  unconventional  definition,  the 
bandwidth  increases  as  the  wall  radius  R  increases  (see  Fig.  6).  If,  however,  the  bandwidth  is 
defined  by  k  range  where  the  growth  rate  falls  off  to  a  prescribed  fraction  of  its  maximum 
value,  as  is  conventional,  it  will  be  much  different  from  that  given  by  Eq.  (69).  This  is  due  to 
the  fact  that  the  maximum  growth  rate  in  the  wavenumber  /c  is  a  decreasing  function  of  the 
wall  radius  R,  as  will  be  shown  later,  thereby  making  the  conventional  bandwidth  a  complicate 
function  of  R. 


B.  Maximixatlon  of  Growth  Rate 


In  this  section  we  will  find  the  dependence  of  the  growth  rate  (A,)  on  the  wavenumber  k 
and  the  wall  radius  R.  Specifically  we  will  examine  the  signs  of  dljdk  and  dA,/d/?  in  -  /? 
space,  and  find  the  position  where  A,  is  absolute  maximum.  This  will  lead  to  the  matching  con¬ 
dition  on  the  wall  radius  R  and  to  the  expression  for  estimated  magnitude  of  the  absolute  max¬ 
imum  growth  rate. 


If  we  differentiate  the  growth  rate  equation  (41)  with  respect  to  k  and  R,  we  obtain 


dk 

BR 


Bk 

BR 


(71) 

(72) 
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where  the  coefficients  and  are  given  by 

i'*- 

III  I  i  \ 

fi'  =  -  (ff?  -  B*')  IS:  -/»•)>  0  i73) 

6>/3  ti 

In  Eq.  (73)  we  made  use  of  B,  >  B  (Eg  (45)),  and  the  guaniiiies  i/and  8^  are  dehned  by 
Egs.  (42)  and  (43).  The  derivaiives,  d<;/dA.  dh/dk.  da/BR,  and  Bh/BR  are  given  in  Egs. 
(53),  (54),  (62),  and  (64).  The  curves  along  which  B\  /Bk  -  0  are  defined  by 

36  j  Bti  jj  n.A\ 

si  “  S  "  ■ 

whose  schematic  solutions  are  illustrated  in  Fig.  8.  Detailed  discussions  on  these  curves'  are 
given  in  Appendix  D.  As  is  shown  in  Fig.  8,  there  are  two  branches  of  B^JBk  -  0  separated 

by  k^  curve.  The  A-dcpendence  of  A,  at  two  representative  values  of  R  are  also  shown  in  Fig. 

8  The  lower  BS  jBk  -  0  curve,  which  gives  higher  maximum  growth  than  that  for  upper 
curve,  intersects  with  the  stability  boundary  k^  and  with  kt  curve  at  the  same  A,^  As  discussed 
in  .Appendix  D,  the  distance  between  the  lower  A-niaximum  curve  and  A,C  line  for  R  <  R„ 
(i.c.,  to  the  left  of  A1  curve  in  Fig.  8)  **  is  given  by 

—  T)**  '  when  €«  ~  »)*.  (75) 

where  cr  =  [Af  -  R„\.  Note  that  6  >  S'^v'e  for  the  stability  curves  Ar,  t*  —  ij'  '  (Eg. 

(70)).  Similarly,  we  can  find  graphic  solution  for  BXJBR  -  0  defined  by 
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In  Fig.  10.  ihe  compleie  k-,  R-  dependence  of  A,  is  shown.  Noie  ihal  dSjBR  -0  curve 
intersect  with  the  lower  dljdk  —  0  curve,  and  its  intersection  coordinates  {R,  /f),  where  the 
growth  rate  is  the  absolute  maximum,  are  given  by  (from  Egs.  (75)  and  (77)) 

T).  (78) 

^  -  t).  (79) 

where  the  expressions  for  R„  and  kH  appear  in  Eqs.  (61)  and  (56)  respectively.  In  view  of 

smallness  of  tj  (Eq.  (50)),  for  practical  purpose,  we  can  conclude  that  the  absolute  maximum 
growth  exists  at 


R  ^  R„- 


y.-  +  I 


P: 

The  corresponding  real  eigenfreguency  <u,  is  then 

w,  = 

The  growth  rate  at  "  R„  and  k  “  can  be  shown  as  (see  Appendix  E) 


A,  - 


>/3 


2v 


y.- 


_ 

y..^M  (y.-  +  iM  s 


(T..)  lix,) 


\/) 


(80) 

(81) 

(82) 

(83) 


where  the  explicite  expression  for  tj  at  R„  is  shown  .  and  the  quantity  X,  is  detined  by 

which  is  Xi  (Eq.  (13))  at  R  -  R„. 


y.+i 


(84) 


In  summary,  we  have  optimized  the  growth  rates  of  the  long  wavelength  modes  with 
respect  to  the  wall  radius  R  and  the  axial  wavenumber  k  for  given  parameters  (/»,  /,  s,  »>). 

The  analysis  has  resulted  in  the  wall  radius  and  the  wavenumber  matching  conditions  (Eqs. 
(80)  and  (81)),  under  which  the  instability  gives  highest  growth  rate,  and  it  has  been  successful 
in  obtaining  approximate  magnitude  of  the  growth  rate.  In  unnormalized  quantities  (Eq.  (10)), 
these  results  are  as  follow. 
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replaced  in  Kgs  (85)  and  (88).  li  is  not  repetitive  to  emphasize  that  the  possibility  of  obtaining 
these  approximate  results  rest  on  the  tenuous  density  (••  «  1.  or  rj  «  1)  assumption. 


V.  PARAMETRIC  OPTIMIZATION  OF  GROWTH  RATE 

In  previous  chapter  we  have  maximized  the  growth  rate  of  the  long  wavelength  mode  with 
respect  to  the  wall  radius  H  and  the  axial  wavenumber  k  while  the  rest  parameters  are  assumed 
to  he  given  In  this  chapter  we  further  optimize  the  growth  rate  with  respect  to  the  beam 
center  location  (section  .A)  the  radial  node  number  n  (section  B).  and  the  azimuthal  mode 
number  /  (section  O  in  that  order.  I'or  each  step  of  optimization  we  start  with  the  growth  rate 
already  optimized  with  respect  to  all  previous  parameters  unless  stated  otherwise.  In  Hnal  sec¬ 
tion  (O)  we  will  discuss  the  dependence  of  the  growth  rate  on  the  magnetic  harmonic  number 
.1.  Again  we  adopt  the  notation  'as  optimum  value. 


A.  Beam  Center  Opilmitatlon  (r„) 

For  given  r,  s,  /  and  n,  the  k  and  R  optimized  growth  rate  is  given  by  Fq.  (83). 

The  only  part  in  Eq.  (83)  that  is  dependent  of  the  beam  center  location  r,  is  3,./(.V„)  term 
through  .Vo  =  .V,,  r„.  Namely, 
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A,  . \J,.,  (89) 

where 

^0  “  *^/«  ^o‘  (8) 

Before  optimization  of  J,-,iX„)  with  respect  to  r„  (Af„),  we  recognize  that  there  exist  a  physical 
constraint  in  variation.  By  insisting  that  the  beam  must  lie  within  the  conductor,  that  is,  in 
unnormalized  real  values, 

C*'  +  '■/'**'  <  /f"*'.  (90) 

which  can  be  expressed  with  normalized  variables  as 

^0  ^  -Kn  ~  Xi_.  (91) 

From  properties  of  Bessel  functions,  we  can  show  that 

1  >  ly._,  >  U.-/(^)|. 

if  s  ^  /.  and  X  ^  (92) 

and 


1  -  ly„(0)|  >  |y.,  (A')|  if  A'  >  0.  (93) 

Here  is  the  first  positive  zero  of  y',.,(.x)  when  s  -  /  ^  0.  It  is  evident  that  the  max¬ 

imum  value  of  (A'„)|  depends  on  (s  -  /),  and  when  s  ^  /  on  the  magnitudes  of  A’,,  and 
We  conclude  that 

Ui-/ (Ar..)L.,  -  (94) 

where 


/•’  = 
A/  — 


1  if  5  -  A 

|y,-;  (A'/,_/n)|  if  S  ?«  /.  Xi„  —  A/  ^  A'|j_;|i. 
|y,-;  (A/,  —  Aft)  I  if  s  ^  /.  Xi„  —  Xi  <  A'|j.;||, 


and  the  corresponding  optimum  beam  center  location  r„  is  given  by 


(95-a) 

(95-b) 

(95-c) 


X,n 


if  s  —  /. 

if  s  ^  /.  Xi„  -  A(  > 


if  S  5^  /.  A';,  —  A';  <  A'|,.)||. 


(96-a) 

(96-b) 


(96-c) 
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Note  (hat 


7;;  (9s-tf)  >  y‘;  (95-/.)  >  j".,  m-c).  (97) 

In  this  regards  the  process  of  r„  optimization  is  in  part  that  of  optimization  of  the  azimuthal 

mode  number  /  for  given  magnetic  harmonic  number  s.  With  r„  optimization  given  by  Egs. 

(94)  to  (96),  the  maximum  growth  rate  A,  becomes. 


A,  - 


2^ 


V; 


•V/ 


A';  A'  (.V/ ) 


1 .1 


(98) 


(y.-  +  1)  I  s 

We  mention  here  that  the  optimized  beam  center  location  (i.e.,  s  -  r,,  “  O'  corresponds  to  so 

called  'rotating  beam',  which  is  quite  different  from  usual  gyrotron  gun,  and  is  best  produced 
by  cusp  magnetic  fields  such  as  Maryland  ERA  machine.  We  also  add  here  that  when  the  beam 
center  optimization  is  not  realized  to  its  potential  as  the  case  is  for  Eq.  (95-c),  there  will  be 
competition  between  the  wall  radius  optimization  (Eq.  (80)),  and  r„  optimization  via  wall  radius 
variation.  This  liind  of  competition  is  present  in  almost  every  step  of  optimization,  for  which 
more  careful  study  is  in  n^ed.  It  is  important  to  point  out  that  for  s  -  /  device,  this  R  -  r„ 
competition  is  absent. 


B.  Radial  Node  Optimization  (n) 


The  dependence  of  the  growth  rate  on  the  radial  node  number  (n),  after  r,.  optimization, 
is  given  by 


A  oc 


1 


I  .< 


(99) 


The  term  originally  appeared  in  the  dispersion  equation  from  orthogonal  integration  of 
Bessel  ’'unction,  that  is. 


,v„  =  j/'"  Xdx  y,-(.v). 


(UX)) 
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Since  Xi„  is  an  increasing  function  of  |/|  and  n,  it  is  obvious  that  Ni„  is  also  an  increasing  func¬ 
tion  of  |/|  and  n.  Therefore  we  conclude  that  the  growth  rate  A/  is  a  decreasing  function  of  the 
radial  node  number  n.  Caution  must  be  taken  here  about  the  competition  between  n  and  r„ 
optimizations.  As  n  decreases,  so  does  X,„,  thereby  raising  possibility  for  violation  of  full 
optimization  of  beam  center  location  (Eq.  (95-c)).  This  competition  is  absent  again  for  s  -  / 
device. 


C.  Azimuthal  Mode  Optimization  (A 

The  dependence  of  the  growth  rate  on  the  azimuthal  mode  number  (A  for  given  magnetic 
harmonic  number  s  is  given  by 


A,  .  W 


1/3 


(101) 


after  the  beam  center  radius  (r„)  optimization  (Eq.  (98)).  Obviously  A,  is  a  function  of  both  5 
and  /,  and  the  combination  makes  it  difficult  to  conclude  on  /  dependence  in  simplicistic  way. 
However,  since  J°,  -  J^-i  (Jr|,_/||)  is  rapidly  decreasing  near  /  -  s,  we  expect  that  A,  will  be 
maximum  when  /  -  s,  unless  the  magnetic  number  s  is  extremely  high.  If  we  ignore  the  full 


optimization  rule  restriction  on  beam  center  radius  (Eq.  (9S-b)),  the  quantity 


‘'In 


1/3 


IS 


completely  independent  of  other  parameters  and  their  numerical  values  are  easily  calculated. 


Fig.  11  gives  the  plot  of  the  quantity 


1  /o2 

nJ- 


1/5 


versus  the  azimuthal  mode  number  /  for 


s  -  1  to  4.  The  values  are  not  to  be  compared  for  different  s,  since  the  quantity  has  not 
included  other  dependences.  However  for  given  s,  these  values  are  proportional  to  their 
growth  rates  (Eq.  (101)),  so  that  the  magnitudes  for  different  /  can  be  directly  compared.  At 
this  low  magnetic  number,  it  is  obvious  that  s  ••  I  growth  rate  is  much  higher  than  that  for 
s  ^  I.  The  points  with  cross  marks  indicate  the  violation  of  the  beam  center  radius  optimiza- 
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That  is.  for  the  azimuthal  wavenumber  /  with  crossed  points,  the  growth  rates  are  less  than 
those  shown  in  Fig.  II.  The  graph  also  reveals  that  except  for  s  -  I.  the  next  highest  growth 
rate  is  for  the  mode  with  /  -  s  -  I.  When  s  -  I.  the  growth  rate  for  /  -  0  is  less  than  that  for 
/  -  2.  The  comparison  of  TE,|,  with  TEm,  (azimuthally  symmetric  /  -  0  model  shows  that 

A,  (r£,„) -2  X  A,  (r£„,.).  (102) 

where  ££,„  notations  are  used.  The  optimization  rule  becomes  then 

h  -  s  (103) 

with 


/j  “  s  —  I  if  j  >  I . 

/j-0  ifj-l, 

as  next  optimum  value.  With  /  —  j  mode,  the  growth  rate  becomes 


(104) 


A 

'  2 


>: 


(>r  +  \)  (  s 
The  corresponding  beam  center  location  is 


J  \  iX,  ) 


I /.I 


(105) 


f„  -  0. 

D.  Harmonic  Number  (s)  Dependence 


(106) 


The  growth  rate,  optimized  on  wavenumber  (A"  =*  A^),  wall  radius  (£  =*  /?„),  beam 
center  at  radius  (r„  -  «).  and  azimuthal  mode  number  (/  -  s)  depends  on  the  magnetic  har¬ 
monic  number  (.v)  according  to  (Eq.  (105)) 


A,  «  1^  ^  (y. 

I  'V|,  .s' 

where  the  argument  of  Bessel  function  is  defined  by 


i/j 


(y, + 


I'j 


s. 


(107) 
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It  is  obvious  from  E«j.  (84)  that  general  properties  of  A,  on  s  are  difficult  to  obtain  It  is,  how¬ 
ever.  possible  to  get  an  approximated  expression  for  )  at  certain  conditions.  As  shown  in 


Appendix  F.  if 


A'  y  ' 

s  7-"  -  «  2. 

(y,  +  1) 


,,  ,1-.^  I  _ Y: 

’  '  2  {2  (y,  -I-  1))'^^ 


where  e  is  the  base  of  natural  logarithms  (c  -  2.71828).  The  condition  (108)  can  be  con¬ 
sidered  as  limits  both  on  the  relativistic  factor  of  beam  electrons  (/Si.  y.)  and  on  the  magnetic 
harmonic  number  s  for  the  validity  of  approximate  expression  (109).  We  mention  here  that, 
for  usual  gyrotrons,  the  relativistic  factor  is  so  low  that  expression  (109)  is  valid  for  all  the 
range  of  magnetic  number  .v  in  interest.  When  the  Eg.  (108)  is  satisfied,  it  is  easy  to  show  with 
(109)  that 


2  fii  yI 


The  other  .^dependent  term  in  A,  (Eg.  (107)),  namely.  1//V„,  is  a  decreasing  function  of  ,s.  but 
as  can  be  seen  from  Fig.  11,  it  is  a  weak  function  of  .v,  compared  with  .v-dependence  of  Fg. 
(110).  We  can  therefore  conclude  that 


A,  (.r)  I 


2(y,  1)  • 


A,  (j)  -  A,  (*  -  n 


. I2(y...,)l 

If  the  relativistic  factor  is  low  such  that 


y!  J_ 

y,  +  1  e' 
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then  the  growth  rate  Er).  (112)  is  a  monotonically  decreasing  function  of  s.  Usual  gyrotron  dev¬ 
ices  easily  satisfy  the  condition  (113).  The  dependence  of  the  growth  rate  of  s  -  /  mode  on  s  is 
illustrated  in  Fig.  12. 


If  we  combine  all  the  optimizing  processes,  we  obtain  the  highest  maximum  growth  rate, 
which  is  given  by 


Z  max  ^ 


fix  y,’ 


yo^w  (y,  +  1)^ 


i/.t 


The  highest  growth  rate  of  Eg.  (114)  is  obtained  with  the  following  optimizations. 


k  =*  A*  -  -j-  (y.  -  1) 

fi; 

vn 


y.-  +  • 


(114) 


(81) 


(80) 


K  -  0  (<)6-a) 

/-w-s-1  (103,111) 

In  unnormalized  guantities,  the  highest  growth  rate  is  for  TE^^\  (/  -  I,  n  -  1,  s  -  I)  mode 

with 


jrenl 


LfCtJ 


y.-  +  ) 


2y/ 

y.--i 


i/j 


c3fii 


fi: 


s<"./y<.  ’ 
i"./y.. 


r.,""  -  0. 

i"*'  -  5  y.  «u,/y„. 

a  real  X^3  j 

A,"*'  -  —  sw,/y„ 


fixf^: 


\li 


(115) 

(116) 

(117) 

(118) 

(119) 


jy..^^M  (y. D’ 

We  note  that  A,""'  (Eg.  (114))  is  a  monotonically  decreasing  function  of  y..  for  given  y„.  That 
is,  for  given  total  particle  energy,  the  growth  rate  is  highest  when  the  electron  has  no  axial 
speed. 


So  far  we  have  considered  the  parametric  optimizations  of  the  long  wavelength  modes. 
This  processes  have  been  eguivalent  to  maximizing  (Eg.  (48)).  We  note  that  the  growth 
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rate  for  the  short  wavelength  (see  chapt.  Ill)  is  proportional  to  and  conclude  that  the 
above  procedures  will  apply  the  same  way  to  the  short  wavelength  modes  as  well  except  for  the 
exponent  (1/2)  instead  of  (1/3)  for  the  long  wavelength  modes.  It  should  be  emphasized  that 
while  the  long  wavelength  mode  possesses  optimum  wavenumber  and  wall  radius,  the  short 
wavelength  growth  rate  increases  as  |A;|  increases  and  as  R  decreases. 

VI.  CONCLUSIONS 

We  have  examined  the  linearized  growth  rates  for  modes  of  the  electron  cyclotron 
maser  instabilities  with  zero-temperature  beam  electrons.  We  have  found  (chapt.  Ill)  that  there 
exist  two  different  types  of  instability  modes  distinguished  by  their  wavelengths;  the  short 
wavelength  mode  and  the  long  wavelength  mode.  Assuming  that  the  electron  beam  is  tenuous, 
we  have  been  able  to  obtain  approximate  expressions  for  the  growth  rates  for  both  types  of 
modes  (Eqs.  (35)  and  (41)),  and  found  that  the  growth  rate  for  the  long  wavelength  mode  is 
much  higher  than  that  for  the  short  wavelength  mode  at  usual  parameter  range.  We  have  also 
found  that  while  the  long  wavelength  mode  has  optimum  axial  wavenumber  and  wall  radius 
that  maximizes  the  growth  rate,  the  short  wavelength  mode  does  not  (see  Fig.  4). 

We  then  have  maximized  the  growth  rate  of  the  long  wavelength  mode  with  respect  to 
various  parameters,  and  have  found  the  optimizing  conditions  on  the  parameters.  These  results 
are  given  in  Table  1. 

The  growth  rate  is  found  to  be  scaled  with  the  harmonic  number  s  accoridng  to  Eq.  (112), 
namely 

A,(j)  -  A,(j  -  1) 
when  the  beam  electrons  are  weakly  relativistic. 


fir  y? 


2  (Y:  +  1) 


5-1 

3 
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We  add  finally  that  some  of  foregoing  optimization  processes  may  compete  with  the  other 
optimizations,  and  a  compromise  between  them  may  result  in.  One  example  of  these  situations 
is  given  in  the  text  following  Eq.  (98).  Whenever  this  competition  rises,  the  original  growth 
rate  equation  (41)  should  be  referred  to.  The  numerical  calculations  of  the  growth  rates  have 
been  carried  out  and  those  results  will  be  reported  elsewhere. 
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Appendix  A.  Roots  near  Local  Extrema 

In  this  appendix,  we  examine  the  roots  of  /(.v)  -  0  near  local  extremum  .v„,  where 
,/''(x„)  -  0.  If  /(.v)  is  Taylor  expanded  near  x,„  then 


/(«)  - ./;.  +  y  fo"  +  -g-  fo"  «’  +  •  •  •. 

where 


In  case  when 


(A-l) 


(A-2) 


»  1*1  Ifo-l 

the  equation  /(c)**0  can  be  approximated  as  a  quadratic  equation.  Namely. 

/(f)  -  ./;.  +  y  f ’  -  0.  (A-4) 

whose  solution  are  given  by 

c  -  ±(  A.  if  >  0.  (A-5a) 

f  -  ±lAl.  if  A’<  0.  (A-5b) 
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I 


Here  is  defined  by 


With  Eq.  (A-5),  the  validity  of  Eq.  (A-4)  (i.e.  condition  (A-3))  can  be  written  as 


\L'V  »  l/J  (A-7) 

In  summarization,  when  xo  is  a  local  extremum  of  function  /fx),  viz  /’(x,,)  -  0.  the  roots  of 

the  equation  /(x)  -  0  are  given  by 


X  -  x„  ±  / 


fo 


f  " 

Jo 


(A-8a) 


if 


0  <  « 

Jo 


f" 


and  by 


X  -  X,.  ±  - 


fo 


2L 


lA-»)a) 


(A-8b) 


if 


/o'" 


« -^  <  0. 

Jo 


(A-9b) 


Notice  that  /(x)  -  0  have  complex  roots  when  and  have  same  signs  (Eq.  9a). 


When  the  condition  (A-7)  is  not  satisfied,  the  quadratic  approximation  Eq.  (A-4)  is  no 
longer  possible,  and  we  have  to  resolve  to  the  cubic  or  further  equation  to  find  the  roots.  The 
condition  (A-7)  can  be  violated  when  |/„"|  is  too  small  and/or  |/„1  is  too  large.  That  is,  in 
order  to  have  complex  roots,  the  equation  fix)  must  have  small  1/„|  large  1/„"1  in  addition  to 
same  sign  of  /„  and 
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Ap^ndlx  B.  Roots  of  CuMc  Equollon 


We  will  summarize  the  properties  of  roots  of  the  cubic  equation  given  by 

x’  +  fix  +  y  -  0.  (B-1) 

where  a,  fi,  and  y  are  real.  With  substitution  of 


.V  =  X  + 


(B-2) 


Eg.  (B-l)  becomes 


where 


i  ’  3/»i’  +  ^  -  0. 


P  =  -  13/J  -  «M. 

9 

9  =  12«’  -  -*■  22yl. 


The  roots  of  Eg.  (B-3)  are  found  to  be 


A*  +  a ' 

*  I  !  J?' 

j  ’  ±  /  ^  a  ' 


where 


i  ~  \  [9  ± 


Here  the  determinant  is  defined  by 


(B-3) 


(B-4) 


(B-5) 


(B-6) 


d' s  +  4p'.  (B-7) 

Note  that  the  necessary  and  sufficient  condition  for  Eg.  (B-I)  to  have  complex  roots  is  given  by 


</’  >  0. 


(B-8) 
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By  making  use  of  Bessel  identity, 

Xl  K(X0  +  -  0. 

Xt. 

the  equation  (C-2)  becomes 

ia. .  ia  [  5^  -  xl  j,(Xl) 

1  Xt  A'(Xt) 

With  substitution  of  Eq.  (C-4),  Eq.  (C-1)  is  now 

^_r  Js(Xt) 

dR~  n,R^  Xt  J‘(Xt)  ■ 

For  derivation  of  Eq.  (64),  we  first  note  that  for  long  wavelength  modes  (i.e. 
k%)  the  term  (Cl}  —  k^),  in  Eq.  (63)  is  positive.  When 

we  have  from  Eq.  (13)  above 


Xt  <  s. 

If  we  denote  the  first  zero  of  A'(X)  for  s  ^  1  as  AT.,  (>  0),  namely.  A'  (;r,,)  -  0, 
Eq.  (C-3) 


-X,] 


A"(X,,) 
AiXsO  • 


The  fact  that  A  has  its  first  maximum  at  -T,,  (i.e.  J',  (T.,)  -  0)  indicates  that 


A(X,0  >  0,  A"  (^,2)  <  0. 

With  Eq.  (C-8),  the  Eq.  (C-7)  becomes 

X,f  >  s. 

By  combining  Eqs.  (C-5)  and  (C-1),  we  obtain 

Xt  <  s  <  X,]. 

Therefore  we  can  easily  show  that 

A  (Xt)  >  0, 

A'  (Xt)  >  0, 
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(C-4) 

(63) 

kt  <  k  < 

(C-5) 

(C-6) 
then  from 

(C-7) 

(C-8) 

(C-9) 

(C-IO) 

(C-11) 
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if  the  condition  (C-5)  is  satisfied.  From  Eqs.  (63),  and  (C-1 1),  it  is  proven  that 

<  0  when  >  y.  (65) 

0  A 

Appendix  D.  dAjdk  ■*  0  and  d^JdR  0  Curves. 


The  curves  along  which  d^Jdk  -  0  are  defined  by 


2  o2 

(74) 

where 

(53) 

(54) 

and 

1  1  1.  1  I  i 

/3^  “  ^  -  Bt)  (Bl  -  Bl)  ^  0. 

(73) 

The  growth  rate  is  given  by 

/T  1-  1- 

A, (fiJ  -Bl). 

(41) 

Along  the  stability  boundary  kg  (i.e.  A,  «  0),  we  have 

/3^  -  0, 

(D-1) 

while  along  a  -  0  curve  (/c* 

in  Fig.  6)  ,  we  find  from  Appendix  E  that 

B\  ~  h.  fl_  =  0. 

(D-2) 

such  that 

-  y 

(D-3) 

On  kH  (i.e.  db/dk  -  0  from  Eq.  (54)),  we  have  from  Eq.  (74) 

a  -  0  or  /3^  -  0, 

(D-4) 
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which  corresponds  to 


k  -  ki  ot  k  ~  kg.  (0-5) 

This  means  (hat  dAJdk  -  0  curve  intersects  with  A*  line  at  kg  and  kt  (see  Fig.  6).  When 

db/dk  we  get  from  Eqs.  (53).  (54)  and  (74), 


which  means  that  the  d^Jdk 
with  difference  of  order  of  tj. 


-  0  curve  is  close  to  o  -  0  (Ar^  in  Fig.  6)  or  to  da/dk 
The  sign  considerations  (Eqs.  (57)-(60))  show  that 


(D-6) 
0  (A  “) 


k  <kl 

dA,  k%  >  k  >  k“ 

Ik  k^>  k  >  k^  ■ 


(D-7) 


\kS  >  k>ki 

which  indicate  that  there  are  two  unconnected  dA,/8A  -  0  curves  separated  by  A"  curve  (Fig. 
7).  For  region  between  the  stability  boundary  (.kg)  and  R„  along  A*  (Fig.  7),  we  have 


da 

dk 


86 

8A 


~  V  «*. 


(D-8) 


where  €g  is  the  /^-distance  of  8A,/8A  -0  curve  from  R„,  and  e*  is  the  A-distance  of 
dAjM  -  0  curve  from  A,f.  If  ««  is  given  by 


than  from  Eq.  (D-8) 


*/» 


(D-9) 


(D-10) 

This  is  interpreted  that  as  8A,/8a  -  0  curve  moves  from  ki  toward  kg,  it  slowly  rises  in  A- 
direciion  and  sharply  decreases  down  to  A*  at  kg  as  shown  in  Fig.  8. 

The  curve  d^,/dR  *  0  is  defined  by 


a 

dR  ^ 


(76) 
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where 


dR 

db  is^-X,')  J,{X,,) 

BR  ~  n,/f^  X,  J.ix,) 

From  Eqs.  (62)  and  ,  (65)  of  Appendix  C,  we  have 


In  particular,  we  note  that 


(D-ll) 


(D-12) 


which  indicates  from  Eg.  (76)  that 


lal-n  (D-13) 

With  sign  considerations  (Eqs.  (57)-(60)).  we  conclude  that  8S,/dk  -  0  curve  lies  to  the  left 

of  k%,  curve  with  difference  ^-direction  in  the  order  of  77.  This  is  illustrated  in  fig.  9. 


Appendix  E.  Growth  Rates  Along  k'‘. 


In  this  Appendix  we  will  obtain  the  expression  for  the  growth  rate  along  A'i  curve  (see 
Fig.  7)  where  a  -  0.  The  growth  rate  A,  is  given  by 

>/3  }  f 

A,-^(fl^-fli).  (41) 


*■  =  I  27  * 

along  li  e.  a  —  0  (Eq.  (15))).  it  can  be  shown  that 


ff*-  6.  ff.-O, 
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and  the  growth  rate  (Eq.  (41))  is  given  by 


Next  we  find  the  expression  for  b  along  .  With  the  substitution  o(  k  -  k%  into  the  expres¬ 
sion  of  b  (Eq.  (52)).  we  obtain  with  the  definitions  (17)  and  (18), 


‘t  2R*  (1  ±  /3,  (1  -  1//?^)'/^) 
Therefore  we  get  the  expression  of  the  growth  rate  on 

A  I  .  -  ^  - a - 

**  2  2/?'' (1  ± /3..  (1  -  l//?^)*''^! 


In  particular,  when  R  —  R„  (Eq.  61) 


y..  +  1 


and  k  -  k,^  kl  (R  -  R„)  (see  Fig.  7),  we  have 


>/3  2t,  y,  1"' 


Here  we  have  made  use  of  the  relation 


1  y..  -  1 


y./3.- 


The  expression  for  t)  is  given  by  Eq.  (48) 


where 


^  ^  U.-/  (XJ  J'  (Art)I^ 

yo”/it 


V  - 


When  R  —  R„  (Eq.  (61)),  the  argument  of  Bessel  function  becomes 


X,  s  s 


--  2/8i^y-^ 


y..  + 1 
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and  7)  (R„)  can  be  now  written  as 


yt(R..)  -  -V  -f  h:  +  1)  IA-,  (X„)  y;  (Xi,)V. 

yo^tn  s 

Upon  substitution  of  Eq.  (E-7),  the  growth  rate  at  /?  •*  R„  (Eq.  (E-5))  becomes 


-t  7  -I  n  —  J^-i^x„)  j:  (kt) 

2  y„N,„  (r.-+l)  s 


Appendix  F.  Asymptotic  Behavior  of 

In  this  appendix,  we  will  find  an  asymptotic  expression  for  Bessel  function  J^(X),  where 
the  argument  X  is  given  in  the  form  of 

X~si.  (C-1) 

Note  that  J^{Xi),  which  appears  frequently  in  the  growth  rate  equation  (e.g.  Eq.  (105)),  is  of 

this  type,  when  Xi^  is  given  by  Eq.  (13)  (or  Eq.  (84)).  The  order  of  Bessel  function  s  is 

assumed  to  be  integer  greater  than  zero.  By  series  representation,  Bessel  function  J^'iX)  can  be 

expressed  as 


j,(x)  j  i  n. 

2  2  ^To 


where  the  /rth  term  T,  is  defined  by 


t-a  -  f  .XU  (2k  +  s)  _X 
’  “  k]  (s  +  *)!  2 


if  we  define  the  ratio  of  consecutive  1 7's*|’s  as  r/,  namely. 


A  -  iT-r' 

r,  =  — : — rr 


;*^'l  (2k  +  S  +  2)  jr 

r*|  “  (*  +  1)  (s  +  *  +  1)  (2k  +  s)  2 


then  it  can  be  shown  that 
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and 


Therefore,  we  conclude  that 


r,  '  >  r, '  if  A|  >  A  j. 


j:{x)  ^  I  ' 


T 

2  (5-1.!  I  2 


if 


(C-5) 


(C-6) 


(C-7> 


(C'-8) 


r;’  «  1 

The  expansion  condition  (C-7)  can  be  rewritten  as,  with  Eg.  (C-1), 


(5  +  1)  ' 

which,  in  turn,  can  be  rewritten,  in  a  stronger  condition,  as 


s  «  4,  (C.9) 

The  condition  (C-9)  may  be  violated  with  large  or  with  large  s.  When  condition  (C-*))  is 

satisfied,  y,'(5()  (Eg.  (C-6)  with  Eg.  (C-D)  becomes 


A  isO  -  j  I  ./; 


(C-IO) 


where  ./,  is  defined  by 


We  note  that 


./;  H  * 


<-i 


(s  -  \)\ 


./;  - 


/i  “  I  if  5  —  I. 

5 


if  5  >  1, 


(C-I2) 

(C-13) 


1(5  -  I)!)'/’- 

and  that  if  5  >  I,  is  a  slowly  increasing  function  of  s,  approaching  e  (-  2,71828)  as  5 
increases  to  infinity.  In  fact  we  can  approimate  /,  as 


./;  =»  (!•)'-' 


(C-I4) 
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for  5  >  I.  The  approximation  (C-14)  is  exact  when  s  •  1,  and  at  s  -  2,  where  equation  (C- 
14)  is  the  worst  approximation  for  /,  (Eq.  (C-lD),  the  difference  is  0.71828  (e-2).  With 
approximate  expression  for  J\  (Eq.  (C-14)),  //(s)  now  becomes 


A'iW  ■  2  1  2  ^ 

We  finally  Id  that  if  Xi  is  given  by  Eq.  (84),  the  expansion  condition  (C-9)  is  now 


(C-15) 


and 


fix  y:  , 


(C-16) 


l(X,) 


1 

2  ‘'l2(y, -i-  1)1'^^ 


(C-17) 
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Tabic  I.  Parametric  Optimization  of  the  Long  Wavelength  Modes. 


Quantity 

axial 

wavenumber 

wall 

radius _ 

beam 


center  r, 


radius _ 

radial 

node  number 
azimuthal 
mode  number 


Optimization  Condition 


Ti-i  iw./y.. 


0. 


R  -  rt. 


if  s  -  In 

if  s  ^  X,„  -  Xi  > 

if  s  ^  I.  X,„  ~  Xl  <  -V|,-,n 


1 


S 
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center  radii  r„  are  contained  in  the  cylindrical  wall  of  radius  R.  The  constant  equilibri* 
um  magnetic  field  is  in  direction  (coming  out  of  paper). 


Fig.  2.  Dispersion  curves  of  the  normalized  waveguide  (cu„)  and  the  Doppler-shifted  cyclotron 
motion  (n,)  versus  the  norlmalized  axial  wavenumber  k.  Also  shown  is  | Ail  lines  for 
reference.  Note  that  the  cutoff  frequency  of  <o,,  is  y./R,  and  that  0,^  -  at  k%.  n} 
■■  o»/i  at  k%,  n,  *  0  at  k„.  The  graph  shown  is  for  the  wall  radius  >  I.  Designated 
as  I  and  II  are  two  wavenumbers  at  which  Figs.  (3a)  and  (3b)  are  drawn. 


cates  the  instability  of  the  long  wavelength  mode 
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Fi»  4.  Schematic  plot  of  the  growth  rate  A,  versus  the  axial  wavenumber  k.  For  k  <  kt  and  k 
>  ikt.  the  short  wavelength  mode  exists,  whose  growth  rate  monontonically  increases 
to  as  I* I  increases  to  inhnity.  For  kt  <  k  <  kt-  the  long  wavelength  mode 

is  present,  whose  growth  rate  has  a  maximum  with  its  magnitude  in  the  order  of 
Note  that  these  two  modes  are  separated  by  two  stable  regions  near  k  -  kt- 
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k-'t,  6b/dk  -  0  at  k^  and  da/dk  -  0  at  k^  (see  Eqs.  (S3)-(60)) 


Fig.  6.  Plots  of  (a^l21)  and  (-b/A)  versus  the  wavenumber  k  for  representative  wall  radii 
designated  from  I  to  V.  Here  Rj  >  Rn  >  R„i  I  >  R^  >  R^.  Instability  exists  in 
/k-space  between  the  two  intersecting  points  of  the  curves  (a'V27)  and  (-6/4)  at  each 
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bility  boundary  on  which  A,  -  0. 
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Fig.  II,  Orowih  rate  normalized  by  U.’t  V;  >1^''’  versus  Ihc  azimuthal  mode  number  /,  for  the 
magnetic  harmonic  number  .r  -  I  -  4.  Only  those  with  the  radial  node  number  «  -  I 
are  shown.  The  values  of  the  growth  rate  are  the  possible  maximum  for  given  (/,  .vl, 
and  for  those  marked  with  open  circles,  the  actual  values  arc  lower  than  those  shown  in 
the  graph  due  to  the  restriction  on  the  beam  center  radius  (see  Ey. 


